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Abstract. We report the first detections of SO2, NS and NO in an extragalactic source, the nucleus of the starburst galaxy 
NGC 253. Five SO2 transitions, three groups of hyperfine components of NO and five of NS were detected. All three species 
show large abundances averaged over the inner 200 pc of NGC 253. With a relative abundance of a few 10 the emission 
of the NO molecule is similar or even larger than that found in Galactic star forming regions. The derived relative molecular 
abundances for each molecule have been compared with those of prototypical Galactic molecular clouds. These results seem to 
confirm that large scale shocks dominate the chemistry of these molecules in the nucleus of NGC 253, ruling out a chemistry 
dominated by PDRs for the bulk of the gas. 
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1. Introduction 

Molecular lines from the central regions of galaxies are among 
the most promising tools to explore and understand the his- 
tory of our universe (e.g. Combes et al. 1 19991 and references 
therein). In many galaxies most of the molecular emission 
stems from highly excited gas within the central few 100 pc. 
However, the processes driving the excitation and the complex 
chemistry of the gas are not entirely clear. 

Molecular studies of nearby active galaxies such as 
NGC 253, M82 or Arp220 suggest that C shocks, photodis- 
sociating radiation, and irradiation by X rays or cosmic rays 
(e.g. Rigopoulou et al. 2002 Giisten et al. 198 1 Farquhar et al. 
I1994> play an important role in the heating and the chemistry 
of nuclear gas. The dominant mechanisms and their relative 
importance are still unclear 

Sulfur-bearing molecules such as SO2 and NS are present 
in a wide variety of interstellar conditions, displaying enhanced 
abundances in the hot cores of high mass star forming regions. 
In hot cores, the chemistry of these molecules is determined 
by grain-mantle evaporation into warm gas (Charnley, 1997 1. 
Abundance ratios such as NS/CS have been suggested to be 
one of the signatures of shocks in hot cores (Viti et al.^001 1. 
The NO molecule is a fairly abundant molecule in the diffuse 
ISM of the Milky Way. This species is the main precursor of 
the N/O chemical network (Half en et al. l2001> . 

In this Letter we report the first detections of SO2, NS and 
NO in an extragalactic source, the nucleus of NGC 253. This 
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region is one of the most prolific sources of molecular emis- 
sion outside the Milky Way (Mauersberger & Henkei n993> . 
The large amount of hot gas (Mauersberger et al. l2003t and the 
high abundance of molecules with special chemistry, such as 
SiO (Garcia-Burillo et al. 2000 1 indicates that the presence of 
large scale shocks should dominate the heating and the chem- 
istry of the molecular gas. Also the very high abundance of 
the three molecules reported in this Letter and their excitation 
can be explained if the chemistry of the nucleus of NGC 253 is 
dominated by large scale shocks. 



2. Observations and results 

Observations of the spectral lines (Table [ij were caiTied out 
with the IRAM 30m telescope. Continuum measurements on 
nearby sources made every ~ 2 hours were used to keep a point- 
ing accuracy better than ~ 3". A wobbling secondary mirror 
was used, with a symmetrical beam throw of 4' in azimuth and 
a switching frequency of 0.5 Hz. 

As spectrometers we used two 256x4MHz filterbanks for 
the 2 mm transitions and two 512x1MHz filterbanks for the 
lines at 1.5 and 3 mm. The SIS receivers were tuned to SSB 
with an image band rejection > 10 dB. Beam sizes were 21" 
(at 3 mm), 17" (at 2 mm) and 10" or 12" (at 1.5 mm). Spectra 
were calibrated with the standard dual load system. Intensities 
are given on a main-beam brightness temperature scale (Tmb)- 
Linear baselines were removed from the spectra and the result- 
ing profiles are shown in Fig.n 
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Fig. 1. Spectra of the SO2, NS and NO transitions observed to- 
wards the nucleus of NGC253 (ajgso = 00H5'"06.0\ 5i950 = 
-25°33'45")- Velocity resolutions are 10, 15 and 20 km s"' re- 
spectively. Ticks indicate the expected position of HF compo- 
nents of NS and NO. The Tmb scale of the y-axis is in mK. 

We detected five transitions of SO2 in the 2 mm band, three 
groups of hyperfine components of NO and five of NS in the 

2, 3 and 1.5 mm atmospheric windows. Profiles are consistent 
with the frequencies in the spectral line catalogues of Lovas 
tl992.) and Pickett et al. (.1998) . We checked that no other 
known molecular transitions, both in the signal and the im- 
age band, significantly contaminate the observed features. The 
only line partially blended is SO2 at 135 GHz that overlaps with 
^'^SO (see Fig. 

3. Analysis 

3.1. Derived line parameters for SO2, NS and NO 

To derive the SO2 line parameters, we have fitted a single 
Gaussian to the observed profiles. Fig. [2 shows the Gaussian 
fits superposed on the observed spectra, and Table ^ presents 
the line parameters derived from the fit. 

Fitting the NS and NO profiles is more complex as their 
ground state presents A-doubling. Thus each rotational level J 
is split into two single rotational levels with opposite parity. 
We use the usual notation to label, in Tabled the lower and 
upper series of single rotational levels as e and / for NS, and 



Table 1. Parameters derived from Gaussian fit to SO2 transi- 
tions and hyperfine structure fitting for NO and NS transitions. 

SO2 



V 


Transition 




^LSR 


Ai'1/2 




(MHz) 




mK km ' 


km s"' 


km s"' 


mK 


134004.8 


8(2,6)-8(l,7) 


380 ( 90) 


242 


91 


3.9 


135696.0 


5(l,5)-4(0,4) 


1000 (200) 


245° 


140° 


5.4 


140306.1 


6(2,4)-6(l,5) 


650 ( 80) 


248 


115 


5.3 


146605.5 


4(2,2)-4(l,3) 


1080 (230) 


241 


148 


6.9 


151378.6 


2(2,0)-2(l,l) 


800 (200) 


225 


158 


4.8 



NS" 

V Transition /TMBtii' Vlsr A»'i/2 Tmb 

(MHz) J~J' mKkms-' kms"' kms"' mK 

115556.2 f-f/ 3400(500) 220° 250° 5.6 

161297.2 2500(700) 199 224 4.3 

161697.2 1 - 1/ 3400(600) 249 278 4.7 

207436.2 |-|« 1500(350) 220° 250° 2.2 

207834.9 5-5/ 1800(350) 220° 250° 2.6 



NO" 



V 

(MHz) 


Transition 
J- J' 


mKkms"' 


Vlsr 
km s"' 


Al'l/2 

km s"' 


Tmb 
mK 


150176.5 
250436.8 
250796.4 


1-in* 
I -In* 
I -In- 


4700 (700) 
3100(320) 
2180 (330) 


244 
250 
215 


157 
150 
147 


13.8 

5.4 

6.1 



" This parameter was fixed in the fit. 

* transitions and integrated intensity refer to the whole group of HF 
components while Tmb and frequency refer to the main component. 



and for NO. Rotational levels are further split into hy- 
perfine (HF) components (Gerin et al. '1992). Fig. ^ indicates 
the location of the HF components, which are unresolved due 
to the broad intrinsic linewidth of the emission from the nu- 
cleus of NGC253. Single Gaussian fits cannot be used in this 
case. We fitted the HF components simultaneously with a comb 
of Gaussian profiles with identical width and the relative fre- 
quency and line intensities fixed to the spectroscopic parame- 
ters of NS and NO. Table[I]shows the derived fine parameters. 

3.2. Column Densities and Rotational Temperatures 

In order to derive reliable physical conditions when compar- 
ing line intensities measured with different beam sizes (Ob), 
we have to make assumptions on the overall extent (Os) of the 
nuclear emission of SO2, NO and NS. According to measure- 
ments of CO and CN transitions with the different beam sizes 
of the 30m telescope and SEST, we derive an extent of the 
emitting region of ~ 20" (240 pc at a distance of 2.5 Mpc; 
Mauersberger et al. 20031 This is in agreement with the size 
obtained by Mauersberger et al. 12003.) from the high angu- 
lar resolution interferometric maps of CS 7 = 2 - 1 of Peng 
et al. ( 119961 . Thus, we can obtain the source averaged bright- 
ness temperatures (T^ - Tmb "^^) over 20" and derive source 
averaged column densities in the upper levels of the observed 
transitions assuming optically thin emission. In the case of NS 
and NO, this was done independently for each series. Fig. |2 
shows the population diagrams for each molecule where the 
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Fig. 2. Population diagrams for the SO2, NS and NO molecules. For NS, the (squares) and 11 series (triangles) are plotted 
and fitted separately. For NO, only the e-series is shown. 



Table 2. Source averaged column densities, rotational temper- Table 3. Fractional abundances compared with similar studies 
ature and abundances. in molecular clouds. 



Molecule N r„, [X]/[H2]'' 
(cm-^) (K) 

SO2 7 10" 14(9) 410"^ 

NS 5 10'3 8(1) 3 10-» 

NO 5 10" 7(2) 3 10^' 

" The uncertainty is a factor of 2 in both directions, due to the uncer- 
tainty of the assumed ^(112)= 1 .7 10^^ cm"^ (Mauersberger et al. . 2003j 



level population can be described by a single rotational tem- 
perature represented as a solid straight line. 

The derived T^ot as well as total column densities and es- 
timated fractional abundances are summarized in Table |2] To 
determine the abundance, we have assumed the H2 column 
density of 1.7 10^^ cm"^ derived by Mauersberger et al. (2003 1 
from their '^CO measurements. The uncertainty of a factor of 
2 in the H2 column density (Mauersberger et al. I2003t does 
not critically affect the conclusions of this work. In general, 
the fractional abundances are high with the most outstanding 
case of NO with an abundance larger than 10"^. For the NS 
molecule, parameters derived for each series are: Nf^s(e) = 
4.410'^cm-2, rrot(e) = 8KandA^Ns(/) = 6.410i3cm-2, 
Ti-otif) - 7.6 K. No substantial difference is found between se- 
ries, in agreement with the result of Gerin et al. (II 9921 1 derived 
from NO in different Galactic molecular clouds. 

4. Discussion 

To establish the mechanism driving the chemistry in NGC 253, 
we compare in Table |3l the measured fractional abundances of 
SO2, NS and NO with those in prototypical molecular Galactic 
clouds dominated by different types of chemistry. We have se- 
lected two dark clouds to illustrate the ion-molecule chemistry 
associated with quiescent gas, three hot cores to illustrate grain 
surface and shock chemistry associated to massive protostars 
and one photodissociation region (PDR) to illustrate the UV 
dominated chemistry produced by OB stars in the main se- 
quence. 

Chemistry dominated by PDRs can be ruled out for 
NGC 253 since the abundances of SO2 and NO are much larger 



Source 


[S02]/[H2] 


[NS]/[H2] 


[NO]/[H2] 




lO-** 


10-' 


io-» 


NGC 253 


0.4 


3 


30 


Dark Clouds 








L134N 


0.4" 


0.2-0.6 * 


20" 


TMCl 


< 0.1° 


0.7-1.2* 


2.7 " 


Hot Cores 








Orion Hot Core 


9.4'' 


0.4" 


30'- 


SgrB2(N) 


3^' 


10^ 


20 ' 


SgrB2(M) 


40' 


0.03 -f 


30 ' 


PDRs 








Orion Bar 


O.OF 




0.2 « 



" Ohishi et al. (T^55); ' McGonagle et al. (1994l; " Gerin et al. 
119921 fT993l : Chamley (19971; ' McGonagle et al. 1 1997 1; ^ source 
averaged Nummelin et al. l2001j : * Jansen et al. 1 1995j . 



than those in Galactic PDRs. It is remarkable that only one 
source in Table |3] the hot core in SgrB2(N), shows larger 
or similar abundances than those observed in NGC 253. This 
is surprising given the expected beam dilution in NGC 253 
for molecules with a high dipole moment like NS and SO2. 
However, it is very unlikely that hot cores dominate the emis- 
sion of the three molecules since the low rotational temper- 
atures derived are inconsistent with the high excitation tem- 
perature (>70K) expected from the dense hot cores. This sug- 
gests that the origin of the large abundance of these molecules 
must be related to a hot core like chemistry, but in regions 
with moderate densities in which molecules are subthermally 
excited. Molecular clouds in the center of the Galaxy show hot 
core like chemistry produced by the ejection of molecules from 
the grains by large scale shocks (see e.g. Martin-Pintado et al. 
200 1'l. Chemistry driven by large scale C shocks has been pro- 
posed to explain the high temperatures and large abundances of 
NH3 and SiO (Garcia-Burillo et al. 120001 Mauersberger et al. 
'2003't. We will now discuss in detail how the abundance of the 
molecules reported in this paper fits into this picture. 
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4.1. SO2 

Towards NGC253, we find a SO2 abundance between that 
found in dark clouds and hot cores. The relatively high rota- 
tional temperature of SO2 rules out that the bulk of its emis- 
sion arises from dark clouds with their typical kinetic temper- 
ature of 10 K. Within the context of hot core chemistry, SO2 
is rapidly formed from the H2S evaporated from grain-mantles 
when r<200K. The same type of chemistry is expected if H2S 
is injected into the gas phase by shocks. The observations of 
SO2 towards SgrB2 by Cummins et al. (1986 1 support this 
idea. They found that SO2 emission arises from two compo- 
nents with very different rotational temperatures, the first one 
associated with the SgrB2 hot core (r,ot=310K) and a second, 
subthermally excited one, with a much lower temperature of 
26 K associated with the envelope. The excitation temperature 
is somewhat larger in the envelope of SgrB2 than in NGC 253. 
However, the latter is averaged over a much larger region. We 
therefore conclude that SO2 arises from the high temperature 
component observed in NH3 by Mauersberger et al. (I2003t . 

4.2. NO 

The large NO fractional abundance averaged over an almost 
200 pc region suggests that NO emission is fairly widespread 
through the whole nucleus of NGC 253. If this emission would 
arise from dark clouds, all the clouds in the nucleus of 
NGC 253 must be of the L134N type, contrasting with what 
we find in the solar vicinity where most of the clouds are 
like TMCl. However, the lack of SiO emission in dark clouds 
(Ziurys et al. 1989 1 suggests that large scale shocks are the ori- 
gin of the NO abundances observed in NGC 253. Like in the 
center of our Galaxy, NO emission from NGC 253 could arise 
from warm moderately dense gas, such as that of the envelope 
of SgrB2 which shows a particular chemistry due to strong 
shocks and/or irradiation by hard X rays (Martm-Pintado T997l 
[2000 1. It is worth noting that models of X ray chemistry (Lepp 
& Dalgarno 1996| predict an enhancement of the NO abun- 
dance (10^^ - 10"^) as observed in SgrB2 and NGC 253. Both 
NGC 253 and the envelope of SgrB2 have enhanced SiO and 
NO abundances suggesting a similar origin. 

4.3. NS 

The NS abundance in NGC 253 is higher than in most Galactic 
dark cloud cores and hot cores. McGonagle & Irvine (1997) 
studied the excitation of NS in several molecular clouds and 
found low excitation temperatures and moderate abundances 
in the hot and very dense regions as expected from the large 
destruction rate of this radical under hot core conditions 
(«(H,) > lO^cm'',!'^ > 100 K). Chemical models of hot cores 
indicate an enhancement of the NS abundance relative to CS 
in the presence of C shocks (Viti et aL lMHl . The high NS/CS 
ratio of 0.4 (with X(CS)=4 lO^^ Mauersberger et al. 2003i 
agrees with the prediction of these models in which a shock 
passed 10 - 20 10^ years after the onset of radiative heating in 
the hot core. This favors large scale shocks as the most likely 
explanation for the NS abundance in NGC 253. It is interesting 



to note that this NS/CS ratio is an order of magnitude above 
the ratio measured by Hatchell et al. (.2002 1 towards six hot 
cores. 

The origin of the large scale shocks is not yet clear. High 
angular resolution images of nuclear SiO emission in NGC 253 
suggest different origins (Garcia-Burillo et al. 2QQQl . In the 
inner circumnuclear disk, shocks from a dense molecular 
outflow driven by massive protostar dominates. In view of the 
low excitation temperatures derived for the new molecules, the 
bulk of the emission is probably not related to this origin. The 
other proposed causes, shocks in the outer Lindblad resonance 
as a consequence of a barred potential or generated by mass 
ejection from the disk, could explain the abundances and 
excitation temperatures. 

Interferometric observations of SO2, NO and NS are 
needed to measure their spatial distribution and to establish the 
origin of the large abundances of these molecules in the nuclear 
region of NGC 253. 
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